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While inflation gives an appealing explanation of observed cosmological data, there are a wide 
range of different inflation models, providing differing predictions for the initial perturbations. Typ- 
ically models are motivated either by fundamental physics considerations or by simplicity. An 
alternative is to generate large numbers of models via a random generation process, such as the flow 
equations approach. The flow equations approach is known to predict a definite structure to the 
observational predictions. In this paper, we first demonstrate a more efficient implementation of the 
flow equations exploiting an analytic solution found by Liddle (2003). We then consider alternative 
stochastic methods of generating large numbers of inflation models, with the aim of testing whether 
the structures generated by the flow equations are robust. We find that while typically there remains 
some concentration of points in the observable plane under the different methods, there is significant 
variation in the predictions amongst the methods considered. 
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I. INTRODUCTION 

The impressive results from the Wilkinson Microwave 
Anisotropy Probe [lj have done much to improve the 
standing of inflation as the leading paradigm for the ori- 
gin of structure in the Universe. However, they have not 
done much in the way of reining in the very large number 
of viable inflationary models, as the uncertainty in the 
key prediction of the spectral index n remains significant 
(and crucially encloses the special value of unity), and 
there is no sign of primordial gravitational waves (speci- 
fied by their ratio r relative to density perturbations). 

The large collection of inflationary models (see Ref. Q 
for extensive reviews) has primarily been developed in an 
ad hoc manner, through selection of potentials motivated 
either by some considerations from fundamental physics 
or by simplicity. Typically these potentials may have sev- 
eral parameters, meaning that the best that observations 
can hope to do is constrain those parameters. Only if 
the potential is particularly tightly defined, for instance 
V((j>) oc 4> 6 , can it be ruled out by present observations. 
Despite the ad hoc way in which the collection of models 
has been constructed, their collective predictions cover a 
fair part of the n-r observational plane, albeit not evenly. 
It does not seem appropriate, however, to interpret this 
as saying that inflation models favour certain types of 
predictions. 

An alternative approach is to throw away the idea of 
taking input from fundamental physics and of enforcing 
simplicity (usually on the potential V(<f>) driving infla- 
tion), and instead seek to generate models of inflation 
via some stochastic process, exploiting numerical tech- 
niques where appropriate. The archetypal such method 
is the inflationary flow equations, introduced by Hoffman 
and Turner [j| and generalized to high order by Kinney 
Q (see also Refs. UISGJ)- Intriguingly, models gener- 
ated via the flow equations exhibit a very clear structure 
in the observational plane, primarily occupying the line 



r ~ or a diagonal locus extending to positive r and neg- 
ative n— 1. The principal aim of this paper is to consider 
whether or not such a structure is a robust prediction of 
stochastically-generated inflation models, or whether it is 
specific to the flow equations implementation. As a by- 
product, we provide a new implementation of the flow 
equations, and also explore the origin of their observa- 
tional prediction more closely. We restrict ourselves to 
single- field inflation throughout. 



II. FLOW EQUATIONS REVISITED 

The flow equations take as their starting point a set of 
differential equations linking a set of slow-roll parameters 
defined from the Hubble parameter H . Following the 
notation of Kinney 3 , these are 
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where primes are derivatives with respect to the 
scalar field. For example, the parameter x Ah equals 
(mpi/47r)ff"/i? and is often denoted r)((j)). Using the 
relation 
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where we define the number of e-foldings N as decreasing 
with increasing time, yields the flow equations 
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where a = 2( 1 Ah) — 4e is a convenient definition. 

As pointed out in Ref. 0, these equations actually 
have limited dynamical input from inflation, since in the 
form d/dcj) they are a set of identities true for any func- 
tion H ((/>), and the reparametrization to d/dN modifies 
only the measure along the trajectories, not the trajec- 
tories themselves. In that light it seems surprising that 
they can say much about inflation at all, but it turns out 
that the flow equations can be viewed as a (rather compli- 
cated) algorithm for generating functions e(</>) which have 
a suitable form to be interpreted as inflationary models 
In the following section we will compare their re- 
sults with more direct ways of generating such functions, 
but in the meantime we will explore the flow equations 
themselves further. 



A. A new numerical implementation of the flow 
equations 

The standard implementation of the flow equations Q 
decides a truncation level for the hierarchy sets ranges for 
the slow-roll parameters within which they are randomly 
selected, and integrates the flow equations either until 
the end of inflation, e = 1, or until a stable late-time 
attractor is reached. In the former case the equations 
are then integrated backwards for a suitable number of 
e-foldings (cither a fixed number such as 50, or one also 
randomly chosen within a range 4j) where the observa- 
tional quantities n, r, and possibly others are evaluated 
and plotted. If a late-time attractor is reached the ob- 
servables are read off at that point. 

Here we use a new and more efficient implementation 
of the flow equations, exploiting the fact that the flow 
equations have an analytic solution discovered in Ref. § . 
This is simply a polynomial in H(ifi): 
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The coefficients Ai can be written in terms of the initial 
values of the slow-roll parameters as 
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We use the ranges specified in Ref. 
those values: 



to randomly choose 



00 



[0,0.8]; 
[-0.5,0.5]; 



2 Ah,o 
3 Ah,o 



-0.05,0.05] ; 
-0.005,0.005] ; 



(7) 



' A 



H.O 



0. 



where the last closes the hierarchy. To allow direct com- 
parison with Kinney's work, we use these ranges through- 
out, though one expects the results to be at least mod- 
estly dependent on the assumptions made here Q . How- 
ever we take the equations to sixth-order, one order 
higher than Kinney's main results, as one of the methods 
we will compare with later can only be implemented for 
even orders. As already shown by Kinney, and separately 
verified by us, such a change in order has negligible im- 
pact on the flow equation predictions. We have carried 
out flow analyses at fifth and eighth orders as well as 
those displayed here. 

Although this solution is analytic, there is still the need 
for one integration in order to determine the number of 
e-foldings as a function of <fi, from Eq. © . However this 
is just a single equation to be integrated, regardless of 
the order to which the flow equations are taken. 

Figure ^ shows the results from the flow equations, 
with the left-hand panel showing our new implemen- 
tation based on the analytic solution, and the right- 
hand panel the traditional multi-equation implementa- 
tion. The same sequence of 40,000 initial conditions was 
used in each case, with most of the points finishing at very 
small r. The values of the spectral index n and tensor-to- 
scalar ratio r were obtained using second-order expres- 
sions and the conventions of Ref. 4]. As expected, the 
diagrams are essentially identical point-by-point, though 
some minor differences occur from the way our imple- 
mentations differ at very small values of e. 

The diagram shows a clear and by now well-known 
structure first noted by Hoffman and Turner Q; in ad- 
dition to the majority of the points at small r, there is a 
swathe of models running in a tightly-defined strip given 
approximately by r = 0.3(1 — n), which is close to but not 
exactly the same as the power- law inflation condition Q . 
While Kinney has been careful not to overinterpret the 
tendency of points to lie in this vicinity, noting that the 
measure on initial conditions is unknown, the results are 
often used to indicate where typical inflationary models 
might lie (e.g. Ref. Q). Our main aim in this paper is to 
investigate the robustness of this structure. 



B. Integration direction 

First however we investigate in a little more detail how 
the structure arises. Kinney's main classification of tra- 
jectories is into those reaching the late-time attractor 
(which all have r — > and n > 1, corresponding to the 
field asymptoting into a non-zero minimum of the poten- 
tial), those trajectories where inflation ends with e = 1, 
and the rejected set of trajectories which are unable to 
sustain sufficient inflation. We make a further division of 
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Flow Equations : Analytic 



Flow Equations : Kinney 



Flow Equations 



FIG. 1: Observables at second order for 40000 initial conditions. The two figures are essentially identical, and reproduce the 
results of Fig. 1 in Ref. Q though at lower point density. 



achieved, in proportion roughly 35 to 1. In Figure [2 
we show only points obtained by forwards integration; 
as there are so few of these we increased the number of 
initial conditions tested up to 600,000, so as to have a 
greater density of points (1225 in all, corresponding to 
just 0.2% of all initial conditions) than in Fig. ^ Wc 
find that those points correspond mainly to the region of 
parameter space where n > I, and do not correspond to 
the main swathe of points seen in Fig.^ That swathe is 
therefore generated entirely from initial conditions that 
have to be integrated backwards in time. The change 
of integration direction is significant, because the critical 
point structure of the system is different in the inverse 
time direction (the flow will typically be to the peaks 
of the potential rather than the troughs). Although the 
backwards integration will pick out those critical points, 
they need not represent a distribution that might have 
arisen from whatever mechanism generated the initial 
conditions. 



FIG. 2: Points in parameter space finishing inflation with 50 
e-foldings or more in the forwards integration. 



those trajectories where inflation ends, into those where 
more than 50 e-foldings of inflation were obtained from 
the initial point (meaning that the location where the 
observables were read off was reached by forwards inte- 
gration from the initial point), and those where less than 
50 e-foldings were obtained from the initial point, so that 
the point corresponding to the observables is effectively 
obtained by integrating backwards in time from the ini- 
tial condition. 

We find that the vast majority of random initial con- 
ditions which end inflation do so before 50 e-foldings are 



However even the backwards integration does not ex- 
plain the flow equations structure, because the only sta- 
ble critical points under backwards integration are at 
r — > and n < 1 . The diagonal swathe is not stable, and, 
as already shown in Refs. 0.[J- the trajectories evolve 
along it, typically following quite complicated paths when 
projected on the n-r plane. Given that the location of 
the swathe is known, it would be nice to explain that be- 
haviour via an analytic argument, but we have not been 
able to find an analytic approximation that makes clear 
why points are able to spend a long time in those parts 
of the plot. 
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III. NEW APPROACHES TO STOCHASTIC 
GENERATION OF INFLATION MODELS 

We now turn to alternative methods of generating ran- 
dom inflation models to contrast with the flow equations. 
As we have seen, the flow equations are equivalent to a 
Taylor expansion of the function H(<fi), which uniquely 
specifies the inflationary trajectory. However it is not the 
only way to do so. Single-field inflation models can also 
be uniquely specified either by giving the potential V(4>) 
or, as noted in Ref. @, by specifying the function e(<j>). 

In order to investigate the robustness of the flow equa- 
tion predictions, one should therefore compare its results 
with those from expansions of these alternative functions, 
as there is no reason to prefer one over another. We 
will consider Taylor expansions of both e(4>) and V ((/>), 
and additionally consider a Pade approximant expansion 
of e(cp). In each case we take the randomly-chosen co- 
efficients to correspond to the same ranges of slow-roll 
parameters used by Kinney for his flow equations analy- 
sis 0. 

If one of the three functions H ((/>), V((j>) and e(</>) is 
specified, then the equivalent form of the others can read- 
ily be obtained. However once the functions are trun- 
cated as expansions at a given order, this correspondence 
breaks down, e.g. a sixth-order truncation of H{4>) does 
not correspond to a sixth-order truncation of e(</>). In 
carrying out these expansions, we are therefore investi- 
gating different subsets of the general slow-roll inflation 
model. Even if the expansions were taken to infinite order 
the correspondence between models will only be obtained 
provided each series is within its radius of convergence, 
which is not guaranteed. On the other hand, by choos- 
ing the same initial values for the slow-roll parameters in 
each case, we are ensuring that the functional forms at 
the initial point are sampled from the same distribution. 



A. e{(f>) as the fundamental input 

The general strategy is similar to the new analytic ap- 
proach to the flow equations, in that we choose coeffi- 
cients randomly to generate a function e(</>), and then nu- 
merically integrate to determine the number of e-foldings 
supported. However at this point we do have to men- 
tion one drawback of using e(</>), which is that when the 
function crosses zero (i.e. into the unphysical region) it 
typically does so linearly which means it does so in a 
finite number of e-foldings. By contrast, the e{4>) gen- 
erated from either the Hubble parameter or the poten- 
tial always approaches zero quadratically, generating an 
infinite number of e-foldings. The e(4>) expansions are 
therefore unable to generate points corresponding to the 
late-time attractor. 



1. Taylor series expansion 
The simplest expansion we can make is a Taylor series 
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where we assume <p = initially. The higher slow-roll pa- 
rameters can all be obtained by differentiating this func- 
tion, and so we can determine the coefficients by ran- 
domly selecting the slow-roll parameters within the same 
ranges as for the flow equations using those relations. 
The relations, taking the expansion to sixth-order, are 
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We are then able to solve the model using a single in- 
tegration to find the iV(</>) relation, as in our new flow 
equations approach. We then apply the same tests to 
the models thus generated: do they lead to an adequate 
number of e-foldings, do they require backwards integra- 
tion to achieve 50 e-foldings, and if satisfactory where do 
they lie in the observational plane? 



5 



2. Pade approximant expansion 

A Pade approximant is an alternative to a Taylor 
expansion, which typically exhibits better convergence 
properties. It is formed of a ratio of two polynomials, 
which may have the same or different orders: 



(10) 



There is a one-to-one correspondence between Pade ap- 
proximants and Taylor expansions of the appropriate or- 
der (K = M + N); those expansions then agree near the 
origin of the expansion but differ as the expansion pa- 
rameter, in this case 0/mpi, becomes of order one, as is 
typical in single-field inflation models. 

We generate the Pade approximants from the Taylor 
series using a routine from Numerical Recipes this 
routine is restricted to N = M which is why we chose 
a sixth-order Taylor expansion above. Having gener- 
ated e(<p) in this way, we proceed as before. One ad- 
ditional caveat is that Pade approximants asymptote to 
constant values; this corresponds to power-law inflation 
but is somewhat artificial and so we exclude points which 
require backwards integration to generate 50 e-foldings, 
and which tend to a constant asymptote between zero 
and one. 



B. V(4>) as the fundamental input 

We can also generate the observables by Taylor ex- 
panding the potential and integrating the e-foldings rela- 
tion to find the values of r and n 50 e-foldings before the 
end of inflation. A similar expansion was used to inves- 
tigate inflationary dynamics in Ref. but that paper 
delineated parameter space regions rather than generat- 
ing ensembles of models. 

We write the potential as 
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The initial condition ranges that Kinney uses are given 
in terms of slow-roll parameters defined from the Hubble 
parameter; to obtain equivalent ranges on the coefficients 
of the potential we use the first-order slow-roll relations 
from Ref. |12|. This gives 
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Those relations are also used to evaluate the observables 
once the position of 50 e-foldings is found. 

C. Results 

In Figure |3 we plot the results for each expansion for 
6000 accepted initial conditions, replotting again the flow 
equations case with that number of points for compari- 
son. The same range has been chosen for the observables 
in each case. 

Each of the expansions is plotted to sixth-order. In or- 
der to check the convergence of the method, we have also 
analyzed the flow equations at fifth-order (to compare 
with Kinney Q) and eighth-order, and the other three 
methods at fourth-order. With the exception of the Pade 
approximant for e(<p), discussed further below, there were 
no significant changes in the distributions obtained indi- 
cating that reasonable convergence had occurred. 

We see that there are significant differences between 
the models, though each does show some level of concen- 
tration in the observable plane. Of the three new meth- 
ods, the Taylor expansion of the potential gives results 
closest to the flow equations, showing indications of the 
same swathe of points with non-zero r, but not however 
reaching to such high values. The classification of points 
is very similar to the flow equations, with 90% trivial 
points, and almost all the remainder requiring backwards 
integration to achieve 50 e-foldings. 

By contrast, the two e(4>) expansions give results which 
are visually quite different. The Taylor series gives a dif- 
fuse ensemble of points, with a preference for n > 1 but 
covering a fairly large fraction of the observable plane. 
The classification of points is also quite different in this 
case, with a higher fraction of points, about 15%, giving 
50 e-foldings of inflation from the forward integration as 
compared to those requiring backwards integration (re- 
call this method does not generate trivial points). 

The Pade approximant expansion of e(4>) gives a dif- 
ferent outcome again, with the separate classifications of 
points leading to different groupings in the plane. The 
vast majority of the points shown correspond to back- 
wards integration. The main feature at to < 1 corre- 
sponds to e(4>) functions which approach zero in the back- 
wards integration, but do generate 50 e-foldings before 
reaching that point; these generate a similar structure as 
the flow equations. The grouping of points to the right 
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Flow Equations 




Taylor expansion for e(0) 



Taylor expansion for V(0) 



FIG. 3: The distribution in the observable plane for each of the four methods discussed in this paper, in each case for 6000 
accepted initial conditions. On the top left are results for the flow equations, to be compared with the Taylor expansion for 
e(4>) (top right), the Pade approximant for e(0) (bottom left) and the Taylor expansion for the potential (bottom right). 



of that at high r corresponds to backwards integration 
points, but this time to functions e(</>) which approach 
one and give 50 e-foldings before that point. There is 
a third distinct grouping, mainly at n > 1 and small 
r, corresponding to points achieving 50 e-foldings in the 
forward integration, but it contains very few points (2% 
of the total). 

However it is less easy to draw firm conclusions from 
the Pade approximant, because we found the method is 
much less well converged than the others. When we ran 
this method at fourth-order rather than sixth-order, the 
same general structures were picked out, but the balance 
of points was quite different with most of the points lying 
in the right-hand set rather than the familiar flow equa- 
tions swathe. By contrast, for the other methods the 



results were essentially unchanged between fourth-order 
and sixth-order. 



IV. CONCLUSIONS 

We have investigated several ways of randomly gener- 
ating sets of inflation models, in order to compare their 
predictions in the observational plane with those of the 
flow equations approach. We have seen that the differ- 
ent methods, all of which are comparably well motivated, 
give significantly different predictions. 

In two of our new methods, we see hints of the struc- 
ture seen in the flow equations, but much less well de- 
fined. Models lying in that region do seem particularly 
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well suited to generating a sufficient number of e-foldings, 
but the narrowness of the band appears to some extent 
to be an artifact of the flow equations implementation. 
In particular, a Taylor expansion of e(<fr), which seems 
as well motivated as the flow equations approach, docs 
not reproduce such a coherent structure in the observable 
plane. 
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